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ABSTRACT 

The oxygen and nitrogen abundance evolutions with redshift of emission- 
line galaxies in the Sloan Digital Sky Survey are considered for four intervals of 
galaxy stellar masses, ranging from 10 11 ' 3 M Q to 10 10 ' 2 M Q . We have measured 
their line fluxes and derived the O and N abundances using recent calibrations. 
The evolution of O and N abundances with redshift clearly shows the galaxy 
downsizing effect, where enrichment (and hence star formation) ceases in high- 
mass galaxies at earlier times and shifts to lower-mass galaxies at later epochs. 
The origin of the scatter in the N/H - O/H diagram has been examined. The 
most massive galaxies, where O and N enrichment and star formation has already 
stopped, occupy a narrow band in the N/H - O/H diagram, defining an upper 
envelope. The less massive galaxies which are still undergoing star formation at 
the current epoch are shifted downwards, towards lower N/H values in the N/H - 
O/H diagram. This downward shift is caused by the time delay between N and O 
enrichment. This time delay together with the different star formation histories 
in galaxies is responsible for the large scatter in the N/H - O/H diagram. 



Subject headings: galaxies: abundances - galaxies: evolution 
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INTRODUCTION 



les. 



Oxygen and nitrogen are key elements in the study of the chemical evolution of galax- 
The N/O - O/H (or N/H - O/H) diagram has been considered by many authors 



Gavilan et al. 



2006 



Edmunds fc Pagelll978tlMatteucci fc Tosilll985l : |Pilvuginlll992l . 11993 : iHenrv et al]|2000 



Izotov et al.ll2006l ). A prominent feature of this diagram is that, in its 



high-metallicity part (12+log(0/H) > 8.3), the N /H abundance rati o shows a large scat- 
ter at a fixed value of the O/H abundance ratio (IHenry et al.l |2000| ; iPilyugin et al.l 120031 : 
Lopez-Sanchez fc Estebanll2010l ). The N/H value for a given O/H contains important infor- 
mation about the heavy element enrichment history of a galaxy, and consequently, about its 
star formation history Therefore, an explanation of the origin of the scatter in the N/H - 
O/H diagram is key to understanding the evolution of galaxies. 

Two main explanations are presently considered to account for the scatter in normal 
spiral galaxies: 1) the time delay between nitrogen and oxygen enrichment; and 2) the 
local enrichment in nitro gen by Wolf-Rayet (WR) stars. Concerning the first explanation, 
Edmunds fc Pagell (119781 ) have noted that, due the fact that O and N are produced in 
stars of different masses, there can be a significant time delay between the release of O, 
produced in high-mass stars, and that of N, produced in intermediate-mass stars, into the 
interstellar medium (ISM). The N/O ratio of a galaxy then becomes an indicator of the time 
that has elapsed since the last episode of star formation. Current stellar evolution models 
predict that N is mainly manufac tured and ejected into the ISM by intermediate-mass stars 
with masses greater than 3-4M^ (IRenzini &: Volil Il981t Ivan den Hoek fc Groenewegenl Il997l ; 



Marigol l200ll ; iGavilan et al. 



2005) 



although massive stars can contribute, at least at low 
metallicities ( IChiappini et al.l 120051 ) . This suggests that the time delay between N and O 



enrichme nt in galaxies is 25 - 400 Myr, depending on the adopted stellar mass-lifetime 
relation ( jRomano et al.ll2005l ). Properties of the observed N/O - O/H diagram can generally 
be reproduced by chemical evolution models of galaxies by taking into account a time delay 



(pVlatteucci fc Tosi 



Gavilan et al. 



19851 : |Pilvugmlll992l Il993l Il999t IHenrv et all 120001 IContini et al.l 12002 



2006). The exact value of the time delay is howe ver still a subject of debate 



(IPilyugin et al.ll2003l : iRicher fc McCaUl 12008 : iThuan et al.ll2010l ). 



Concerning the second explanation, it is known ( ILopez-Sanchez fc Estebanll2010l and 



references therein) that some galaxies with WR features have a hig h N/O ratio. This s uggests 
that the ejecta of WR stars may locally enrich the ISM in N. IHenrv et al.l (120001 ) noted 
that most points in the N/O - O/H diagram cluster at relatively low N/O values. This 
led them to conclude that the lower envelope is the "equilibrium" unperturbed locus and 
that the observed scatter is the result of intermittent i ncreases in N. caus ed by the local 
contamination by WR stars or luminous blue variables. Ilzotov et al.l (120061 ) also arrived at 
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the same conclusion. 



Selective heavy element loss through enriched galactic winds may also introduce scatter 
in the N/O ratios. It is believed, however, that galactic winds play an important role in the 
chemical evolution of only dwarf galaxies, not of giant spirals (e.g. iPilyugin et al.l 120041 ) . 



In recent years, the number of good-quality spectra of emission-line galaxies has in- 
creased dramatically due to the completion of several large spe ctral surveys, and in partic- 
ular of the Sloan Digital Sky Survey (SDSS) (jYork et al.l 120001 ) . This opens the possibility 
of using the SDSS spectral base t o study the evolution of O and N abundanc es in galaxies 



in th e redshift range z < 0.4 (e.g. iTremonti et al.l 120041 ) . In a previous study ( IThuan et al. 



20101 ). we have considered the evolution of O and N abundances in galaxies with different 
stellar masses. We have found clear evidence for galaxy downsizing, where the sites of active 
star formation and hence of metal enrichment shift from h i gh-m ass galaxies at early cosmic 
times to lower-mass systems at later epochs (jCowie et al.l Il996j ) . All enrichment ceases in 
the most massive galaxies at late cosmic times. This downsizing effect provides a remark- 
able opportunity for clarifying the origin of the scatter in the N/H - O/H diagram. Indeed, 
if the scatter is caused by a time delay between O and N enrichment in galaxies, then it 
should be minimized for the most massive galaxies where enrichment has stopped. On the 
other hand, one should expect that the N/H ratios of the less massive galaxies (those with 
significant star formation at the current epoch) to be shifted towards lower values relative 
to the massive galaxies in the N/H - O/H diagram, because there is a time delay between 
N and O enrichment and N has not yet been released in the ISM. 

We examine in this paper whether the time delay is in fact responsible for the large 
scatter in the N/H - O/H diagram. Our previous study was based on the MPA/JHU 



catalo g s of automatic l i ne flu x measurements of the SDSS spectra (see iBrinchmann et al. 



( 120041 ); ITremonti et al.l ( 120041 ) and other publications of those authors). The accuracy of 
such automatic line flux measurements seems is not good enough for investigating the origin 
of scatter in the N/H - O/H diagram. We have thus decided to manually measure the line 
fluxes of the SDSS spectra. This provides us with more accurate line flux measurements, 
especially for spectra of objects with the largest redshifts, those with redshifts between 0.3 
and 0.4. 



2. SAMPLE SELECTION 



We first extract from the MPA/JHU catalogs four subsamples of emission-line galaxies, 
each covering a different interval of galaxy stellar mass, centered respectively on the mass 
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values M s = 1O 1L3 M , 1O 11O M , 1O 1O ' 6 M , and lO ia2 M . To have a reasonable number of 
galaxies in each subsample and in order for the galaxies to cover the whole redshift range in 
a roughly uniform fashion, the mass interval dM was chosen to be small and to vary from 
0.0015 to 0.10 dex, depending on the redshift and Ms- For each mass interval, we have 
extracted from the MPA/JHU catalogs emission-line galaxies with automatic measurements 
of fluxes in the H/3, Ha, [O ll]AA3727,3729, [Om]A5007, [Nll]A6584 emission lines. We 
impose the additional restrictions that the galaxies have equivalent widths EW(H(3) > 10 
A (this criterion applies only to galaxies with redshifts < 0.3, as we did not wish to reduce 
further the very small number of galaxies with z > 0.3), EW([0 ll]AA3727,3729) > 3 A, 
£W([Om]A5007) > 3 A, and £W([Nn]A6584) > 2 A. These additional restrictions insure 
that the chosen galaxies have a reasonably high star formation rate, giving rise to reasonably 
strong emission lines that can be measured with good accuracy. 

Each so chosen spectrum was then examined visually and the noisy spectra were re- 
jected. The line fluxes of the remaining galaxies were then measured with IRAF B. The 



[Om]A5007)/H/3 vs [Nn]A6584)/Ha diagram was then used to reject AGNs ( IBaldwin et al. 



19811), with the dividing line between H II regions ionized by star clusters and AGNs taken 



from iKauffmann et al.l (120031 ). 

The wavelength range of the SDSS spectra is 3800 A - 9300 A, so that for nearby 
galaxies with redshift z < 0.023, the [O ll]A3727+A3729 emission line is out of the observed 
range. Thus, all galaxies in our total sample have redshifts greater than ~ 0.023. For distant 
galaxies with redshift z > 0.33, the [SlljAA 6717,6731 emission lines are out of the observed 
range. The sulfur line intensities are usually used as indicators of the electron density. They 
serve also to distinguish between h ot and warm Hll regions, in the framework of the ON 



calibrations of iPilyugin et al.l (j2010j). Since the sulfur lines are not present, we will assume 
the electron density to be equal to 100 cm -3 and that our sample does not contain hot H II 
regions (i.e. those with 12+log(0/H) < 8.0). The redshift z and stellar mass Ms of each 
galaxy were taken from the MPA/JHU catalogs. When the stellar mass of a galaxy with 
redshift z > 0.3 is not available in the MPA/JHU catalogs, we estimated it from its absolute 



1 IRAF is distributed by National Optical Astronomical Observatories, which are operated by the Associ- 
ation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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magnitude M Z} its colour m z -m r and its H/3 equivalent width, using the relation 



log M s = -0.459 M z - 0.775 (m z - m r ) + 0.363 + 0.083 log EW{Rf3) 
for (m z — rn T ) < —0.3 



(1) 



log M s = -0.61 M z + 0.43 [m z - m r 
for (m z — m r ) > —0.3 



2.58 - 0.24 log EW (H/3) 



where the galaxy stellar mass Ms is in units of solar masses. This relation was derived 
by fitting the data for galaxies in the redshift interval 0.3 < z < 0.4 with available stellar 
masses in the MPA/JHU catalogs. We caution that this relation is not a general relation 
for the determination of galaxy stellar masses. It applies only within the small particular 
redshift interval defined above. We stress also that the derived masses serve only to select 
galaxies with various masses to populate our different subsamples and are not used in any 
calculation. 

Our final database consists of 221 spectra of galaxies with Ms = 10 1L3 M© (with 74 
objects with Ms derived from Eq. 1), 259 spectra of galaxies with Ms = 1O 1LO M (60 
objects with Ms from Eq. 1), 244 spectra of galaxies with Ms = 1O 1O ' 6 M (32 objects 
with Ms from Eq. 1), and 152 spectra of galaxies with Ms = 10 ia2 M Q (7 objects with 
Ms from Eq. 1). The measured emission-line fluxes are then corrected for interstellar 



reddening using the theoretical Ha to H/3 ratio, in the same way as in iThuan et al.l ( 120101 ) . 
O and N abundances a re then estimated for each galaxy, using the recent ON calibrations 
of iPilyugin et al.l ( 120101 ) . These calibrations giv e both O and N abun dances over the whole 
metallicity range with a satisfactory precision. IPilyugin et al.l (120101 ) found that the mean 
differences between O abundances determined from these calibrations and from the direct T e 
method based on the [Om]A4363 line is only ~0.075 dex for O abundances, and only ~0.05 
dex for N abundances. 



3. GALAXY DOWNSIZING AND THE ORIGIN OF THE SCATTER IN 

THE N/H - O/H DIAGRAM 

We first check the accuracy of our abundance determinations. FigJT] shows the N/H 
- O/H diagram. The gray points represent galaxies in our four SDSS subsamples with O 
and N abundanc es derived from the ON calibrations. The dark triangles show H II regions in 
nearby galaxies ( IPilyugin et al.ll2010l ). with O and N abundances derived with the T e method, 
commonly thought to be the most accurate one. It is seen that the galaxies in our SDSS 
subsamples occupy the same region in the N/H - O/H diagram as the Hll regions in nearby 
galaxies with O and N abundances derived through the T e method. This confirms that our 
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ON-calibration-based O and N abundances are quite reliable. The remarkable feature to 
note in the N/H - O/H diagram is the large scatter in N/H values at a given O/H value. 
This holds as much for SDSS objects as for H II regions in nearby galaxies. 

We next investigate the redshift evolution of abundances in each of the four galaxy 
mass ranges (Figj2]). In each panel labeled by the galaxies' masses, the filled gray circles show 
individual galaxies. The solid line is the best least-squares fit to those data. Inspection of 
the upper panel shows that there is no systematic variation with redshift of the O abundance 
in galaxies with masses ~ 1O 11,3 M up to z = 0.4. This implies that these galaxies have 
reached a high astration level some 4 Gyr ago, and have been somewhat "lazy" in their 
evolution afterwards. The average value of the O abundances for galaxies with stellar mass 
Ms = 1O 1L3 M is remarkably close to t he O abundance in the Orion nebula (12+log(0/H) = 



8.51) obtained by lEsteban et al.l (12004 ). I t is also in good agreeme nt with the O abundances 



in nearby luminou s gala xies obtained by iMoustakas et al.l (120101 ). using the calibration of 



Pilvugin & Thuanl (120051 ) . 



Comparison between the different panels of FigfS] shows clearly the effect of galaxy 
downsizing. Massive galaxies do not show O enrichment: they already reach a relatively 
high O abundance at z = 0.4, and that abundance remains nearly constant until z = 0. The 
less massive galaxies do show O enrichment: they have lower O abundances at z = 0.4, but 
these increase from z = 0.4 to z — 0, and at the present epoch, they are nearly the same as 
in massive galaxies. 

We next investigate the origin of the scatter in the N/H - O/H diagram. The upper left 
panel of Figj3] shows that diagram for the subsample that contains our most massive galaxies, 
those with masses ~ 10 1L3 M Q . It is seen that the scatter in the N/H - O/H diagram for 
this subsample of galaxies is small, especially for galaxies with z < 0.3 (see Fig. 4, to be 
discussed later). Using this subsample, we have derived the following relation between N/H 
and O/H abundances 

log(N/H) = 2.596 (±0.118) log(0/H) - 14.359 (±1.001). (2) 

This relation has been derived in an iterative manner. First, a least-squares fit was obtained 
for all data points with z < 0.3. Then, objects with deviations larger than 2a were rejected 
and a new least-squares fit derived. The final fit is obtained when, after several iterations, 
the fits for two consecutive ones coincide. The scatter in the N/H values relative to the final 
fit is cr = 0.037 dex. The derived N/H - O/H relation is shown in the upper left panel of 
Figj3] by a solid line, while the dashed lines show ±2a shifts from that relation. It should 
be noted that this N/H - O/H relation is derived for a relatively small range of O/H and, 
consequently, its slope cannot be determined with a high precision. Therefore, it should not 
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be used for 0/H values beyond the above specified range. To establish a more general N/H - 
O/H relation, a sample of objects covering a more extended range in O/H would be needed. 

The above relation links N/H and O/H in galaxies where O and N enrichment has 
already ceased. One would expect that galaxies with significant star formation in the near 
past would be displaced towards lower N/H values in the N/H - O/H diagram relative to 
this relation, because of the time delay between N and O enrichment. Comparison between 
the four panels of FigfJ] clearly shows this effect. The most massive galaxies define an upper 
envelope in the N/H - O/H diagram and they do not show significant star formation in the 
last 4 Gyr (FigfSJ). The lower- mass galaxies are shifted towards lower N/H values and they 
show significant star formation during the same period (Figj2]). This is solid evidence that 
the differences in N/H for galaxies with a given O/H are caused by the time delay between 
N and O enrichment. 

FigJH shows how the deviations of the N abundance from the N/H - O/H relation 
(Equation 2) vary with mass and redshift. The gray points represent individual galaxies. 
The dashed lines show ±2<r deviations from that relation. The N evolution with redshift in 
galaxies of different masses is a clearly due to a galaxy downsizing effect. As the universe 
ages, the sites of star formation shift from high-mass galaxies (log M = 11.3), where star 
formation and O and N enrichment have ceased more than 4 Gyr ago, to lower-mass galaxies 
(log M < 11.3) where active star formation and O and N enrichment are still ongoing. 
Comparison of FigsH] and [2] shows that the downsizing effect is more prominent for N than 
for O. This is because N is a secondary element while O is a primary one. 

To have comparable numbers of galaxies with low redshifts (0.02 < z < 0.15) in the 
highest-mass subsample of galaxies (Ms = 10 u ' 3 Mq) and in the lowest-mass subsample (Ms 
= 10 10,2 M Q ), we had to use a mass interval dM = 0.10 in the first case, and a value about 
66 times lower, dM = 0.0015, in the second case. This implies that, at the current epoch, 
star formation events occur nearly two orders of magnitude more often in low-mass galaxies 
than in high-mass ones. Because of the time delay between O and N enrichment, there is a 
downward shift towards lower N/H values for lo w-mass galax i es in t he N/H - O/H diagram 



(FigHJ). This naturally explains the finding of iHenry et al.l ( 120001 ) that most of the local 



dwarf star-forming galaxies cluster at relatively low N/O values in the N/O - O/H diagram. 



4. CONCLUSIONS 



The O and N abundance evolutions with redshift in emission-line galaxies from the 
Sloan Digital Sky Survey (SDSS) are investigated. To investigate the galaxy downsizing 
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effect, we have measured line ffuxes of selected galaxies in four intervals of galaxy stellar 
masses, ranging from 10 11 ' 3 M & to 10 10,2 M m . The O and N abundances have been derived 



using the recent accurate calibrations of iPilyugin et al.l (120 IQh . 



We have found that the O and N abundance evol utions with redshi ft of the galaxies in 
our sample clearly show the galaxy downsizing effect (ICowie et al.lll996l ). where enrichment 
(and hence star formation) is shifted from high-mass galaxies at earlier cosmic times to 
lower-mass galaxies at later epochs. 

We have examined the origin of the scatter in the N/H - O/H diagram. We have found 
that the most massive galaxies where O and N enrichment and star formation has already 
stopped, occupy a narrow band in the N/H - O/H diagram, defining an upper envelope. On 
the other hand, the less massive galaxies with significant star formation at the current epoch 
are shifted downwards, towards lower N/H values. The scatter of N/H for a given O/H 
is caused by the time delay between N and O enrichment and the different star formation 
histories in different galaxies. 
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12+log (O/H) 



Fig. 1. — The N/H - O/H diagram. The gray (light-blue in the electronic version) points 
represent individual galaxies from our four SP SS subsamples with oxygen and nitrogen 
abundances derived from the ON calibrations of iPilyugin et al.l ( 120101 ) . The dark (black in 
the electronic v ersion ) triangles show H II regions in nearby galaxies (from a compilation of 



Pilyugin et al.l (120101 ) with oxygen and nitrogen abundances derived with the T e method. (A 



color version of this figure is available in the online journal) 
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Fig. 2. — Oxygen abundances as a function of redshift for the four subsamples of galaxies 
with different stellar masses. In each panel, labeled by the galaxy stellar mass, the gray 
(light-blue in the electronic version) points represent individual galaxies. The solid (black in 
the electronic version) line is the least-squares best fit to those data. The dark (red in the 
electronic version) filled triangle in the top panel shows the Orion Nebula. (A color version 
of this figure is available in the online journal) 
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Fig. 3. — N/H - 0/H diagrams for the four subsamples of galaxies with different stellar 
masses. In each panel, labeled by the galaxy stellar mass, the gray (light-blue in the electronic 
version) points represent individual galaxies. The solid (black in the electronic version) line 
in the left upper panel is the derived N/H - O/H relation (equation 2) for the most massive 
galaxy subsample, with masses ~ 1O 1L3 M . The dashed (black in the electronic version) 
lines show the ±2<r shifts from this relation. These dashed lines are reproduced in all other 
panels. (A color version of this figure is available in the online journal) 
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Fig. 4. — Deviations of the nitrogen abundance from the N/H-O/H relation, derived for the 
most massive galaxy subsample (equation 2), as a function of redshift for all four subsamples 
of galaxies with different stellar masses. In each panel, labeled by the galaxy stellar mass, the 
gray (light-blue in the electronic version) points represent individual galaxies. The dashed 
(black in the electronic version) lines show the ±2cr deviations from the N/H-O/H relation. 
(A color version of this figure is available in the online journal) 



